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We present a comprehensive study of r-process element nucleosynthesis in the ejecta of com¬ 
pact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of 
the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including 
viscosity, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment clo¬ 
sure neutrino transport scheme. The investigated cases are guided by relativistic double neutron 
star (NS-NS) and NS-BH merger models. Our nucleosynthesis analysis includes the dynami¬ 
cal ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the 
relic BH-torus systems. While typically ~ 20 — 25 % of the initial torus mass are lost by viscously 
driven outflows, neutrino-powered winds contribute at most another ~ 1 %. Since BH-torus ejecta 
possess a wide distribution of electron fractions and entropies, they produce heavy elements from 
A 80 up to the actinides, with relative contributions of A > 130 nuclei being subdominant. The 
combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly 
well for A > 90. Varying contributions of the torus ejecta might account for observed variations 
of lighter elements with 40 < Z < 56 relative to heavier ones, and a considerable reduction of 
the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might 
explain the composition of heavy-element deficient stars. 
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Figure 1: Schematical overview of our study. See main text for details. 

1. Motivation and Overview 


The astrophysical site(s) for the production of r-process elements is (are) still mysterious. Re¬ 
cent studies suggest that instead of or in addition to the long favored core-collapse supernovae 
the outflows from compact binary (i.e. double neutron-star, NS-NS, or neutron-star black-hole, 
NS-BH) mergers could provide suitable conditions to allow for the r-process. Most existing nucle¬ 
osynthesis studies [1-6] only consider the prompt (dynamical) ejecta launched during the merger 
phase in a NS-NS configuration. However, also the dynamical ejecta from NS-BH mergers as well 
as the ejecta stemming from the remnants of both merger types need to be taken into account in 
order to understand the full chemogalactic fingerprint of compact object mergers. In both merger 
types, the merger remnant can consist of a BH surrounded by a more or less massive torus of 
neutron-star debris. In [7] we conducted a comprehensive analysis of the nucleosynthetic output 
associated with the ejecta from the merging phase of NS-NS and NS-BH binaries and with the 
outflows of the subsequent long-time evolution of relic BH-torus systems. We summarize the in¬ 
dividual steps of this study in Fig. 1 . The compact binary mergers are simulated with a relativistic 
smooth-particle-hydrodynamics code [8], while the BH-torus modeling is performed with a Eule- 
rian finite-volume Godunov-type scheme [9], supplemented by a shear-viscosity treatment with a 
Shakura-Sunyaev a-prescription for the dynamic viscosity [10]. For the long-time evolution we 
employ a pseudo-Newtonian approximation of the gravity potential for the rotating relic BHs [11]. 
In the time-dependent BH-torus modeling we apply, for the first time, detailed energy-dependent 
and velocity-dependent 2D neutrino transport based on a new two-moment closure scheme [12], 
which allows us to determine the neutrino-driven wind and the neutron-to-proton ratio in the disk 
outflows with higher accuracy than in previous simulations. All simulations include microphysical 
treatments of the gas equation of state (EOS). 

Our nucleosynthesis calculations are carried out in a post-processing step of the ejecta pro¬ 
duced by the hydrodynamical models, using a full r-process network including all relevant nuclear 
reactions [13-15]. Eor the combined analysis we pick cases from larger sets of NS-NS merger 
models, which lead to prompt or slightly delayed BH formation, and NS-BH merger models on the 
one side and BH-torus models on the other side such that the macroscopic system parameters (BH 
and torus masses and BH spins of the merger remnants) match roughly on both sides. 
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Figure 2: Left: Rest-mass density in the equatorial plane during the merger of a 1.4Mq NS with a 5.1 M© 
BH. Dots denote particles (projected into the equatorial plane) which are or will become gravitationally 
unbound while arrows indicate the corresponding velocities. The visualization tool SPLASH was used to 
convert SPH data to grid data [16]. Right: Abundance distributions as functions of the atomic mass for 
the dynamical ejecta of three NS-BH merger cases. Each binary system is characterized in the legend by 
the EOS used in the simulation and the mass (in Mq) of the NS and BH, respectively. All distributions are 
normalized to the same A = 196 abundance. The dotted circles show the solar r-abundance distribution [17]. 

2. Ejecta from the Merger Phase 

Consistent with previously published relativistic results [5, 18, 6], the “dynamical” ejecta of 
NS-NS and NS-BH mergers, which are expelled within milliseconds of the collision of the two 
binary components, were found to possess similar average properties, namely expansion velocities 
of 0.2-0.4c, electron fractions below ~ 0.1, and entropies per baryon of a few k^. The considered 
NS-NS mergers produce ~0.004-0.021 Mq of ejecta, whereas the NS-BH mergers eject signifi¬ 
cantly larger masses, 0.035-0.08 Mq, with very low entropies (;$ I^b per nucleon), because this 
matter is not shock heated as in NS-NS collisions, but originates mostly from the outer tail of the 
tidally stretched NS at its final approach to the BH. Mass lost in NS-BH mergers is also expelled 
much more asymmetrically than in the case of NS-NS mergers (see left plot in Fig. 2): Correspond¬ 
ing hemispheric asymmetry parameters (mass difference between dominant ejecta hemisphere and 
opposite hemisphere, divided by total ejecta mass) are a few per cent for symmetric NS-NS mergers 
and 15-30% for strongly asymmetric ones, but 0.93-0.98 for NS-BH mergers. 

Since the high neutron excess, thermodynamic properties, and expansion timescale are very 
similar, the ejecta of NS-NS mergers as well as those of NS-BH mergers are sites of robust produc¬ 
tion of r-nuclei with A ^ 140 and abundances close to the solar distribution (see right plot in Fig. 2 
for typical abundance pattern resulting in the NS-BH case). This result holds basically indepen¬ 
dently of the considered nuclear EOS and the exact binary parameters and confirms fhe findings of 
previous sfudies based on relafivisfic merger simulafions [5, 19, 6]. 

3. Ejecta from the BH-Torus Phase 

The relic BH-forus sysfems lose mass in neufrino-driven baryonic winds [20] and in oulflows 
friggered by viscous energy dissipafion and angular momenfum fransporf [21-23]. We find fhaf fhe 
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Figure 3: Left: Snapshot of a BH-torus model with a 3Mq central BH with spin parameter 0.8 and an 
initial torus mass of 0.3Mq at about 50 ms after the merger. On the left side the net neutrino heating rate 
due to j3-processes and neutrino pair annihilation is color coded and overlaid with arrows indicating the 
energy-integrated flux density of electron neutrinos multiplied with the squared radius. On the right side the 
density is color coded and overlaid with velocity arrows. Right: Abundance distributions as functions of the 
atomic mass for three BH-torus systems with BH masses of Mbh = 3,4 and 6 Mq and the same 0.3 Mq tori. 
All distributions are normalized to the same solar A = 130 abundance. The dotted circles show the solar 
r-abundance distribution [17]. 


viscously driven ejecta amount up to, fairly model independent, 19-26% of the initial torus mass 
(for our BH spins of Abh = 0.8) and dominate the neutrino-driven ejecta by far. The neutrino-driven 
wind depends extremely sensitively on the BH mass and on the torus mass, which determines the 
neutrino luminosities. The neutrino energy loss rates can reach several 10^^ ergs^' up to more 
than 10^^ erg s^' for each of Vg and Vg for a few 100ms. The neutrino-driven ejecta carry away 
up to about one per cent of the initial torus mass, but their mass can also be orders of magnitude 
lower. The maximum masses that can be associated with ejecta driven by neutrino heating are 
2.5-3.5 X10^^ Mq in the case of high neutrino luminosities for typical durations of fractions of a 
second up to ~1 s. This number agrees with the lower limit of the expelled mass computed by 

[24] for a hypermassive NS as merger remnant, and the corresponding mass-loss rate is roughly 
compatible with the values obtained for neutrino winds of proto-neutron stars in supernova cores 

[25] . Neutrino heating, however, inflates the outer layers of the torus and has a positive feedback 
on the viscously driven mass ejection on the level of several per cent of the torus mass (or up to 
~20% of the viscous-outflow mass). 

Neutrino-driven torus winds exhibit characteristic properties which distinguish them from the 
viscously triggered outflow, namely the tendency of higher mean entropies, higher electron frac¬ 
tions, and larger expansion velocities. Moreover, the neutrino wind is strongest at early times and at 
intermediate latitudes (around 45° away from the equatorial plane), and the entropy, electron frac¬ 
tion, and velocity exhibit a strong pole-to-equator variation with higher values towards the poles. 
In contrast, viscously driven ejecta develop on longer timescales, are more spherical, and their 
properties vary little with angular direction. 
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Figure 4: Left: Abundance distributions as functions of the atomic mass for three combined systems (merger 
model plus remnant model) corresponding to models with the indicated torus masses. All distributions are 
normalized to the same solar A = 196 abundance. The dotted circles show the solar r-abundance distribution 
[17]. Right: Comparison between the HD88609 and HD122563 elemental abundances (in log£ scale) and 
those estimated for the combination of a merger model and a BH-torus model, when the ejected mass of 
the merger model is assumed to be 100 times smaller than the one coming from the BH-torus model. The 
calculated distribution is normalized to the Sr abundance of HD88609. 


Because of their greater neutron excess, viscously-driven ejecta allow for a much stronger 
r-process than the neutrino wind. The combination of both components is far dominated by the 
viscous contribution and matches the solar abundance pattern for all nuclear mass numbers A ^ 90 
fairly well. The abundance pattern for A < 132 is comparatively uniform, but the strength of the 
third abundance peak decreases with higher BH mass (see Fig. 3), and the relative yields of low- 
mass (A 130) and high-mass (A ^ 130) components depend on the value of the dynamic viscosity 
and the detailed treatment of the viscosity terms in the hydrodynamics equations. 


4. Combined Ejecta and Observational Implications 

The consistently mass-weighted combination of the prompt ejecta from the merger models and 
the secular ejecta from the BH-torus models can reproduce amazingly well the solar r-abundance 
pattern in the range 90 ^ A < 240 and therefore also that seen in ultra-metal-poor stars (cf. left plot 
in Fig. 4). In particular, the BH-torus outflows are able to well fill the region A 140, where the 
prompt merger ejecta underproduce the nuclei (but see [6]). Since the relative yields of the relic 
BH-torus systems for A < 130 and A > 130 nuclei depend sensitively on the system parameters, 
whereas the A ^ 140 species are created with a robust pattern during the binary merging phase, we 
expect a larger variability in the low-A regime than for high mass numbers. 

Another interesting possibility is connected to the fact that the mass ejection during NS-BH 
mergers shows extreme spatial asymmetry (corresponding to asymmetry parameters ^ 0.95) but 
the mass loss from their BH-torus remnants is much more isotropic. This can lead to a strong 
suppression of the dynamical ejecta component relative to the torus outflow in observer directions 
pointing away from the hemisphere that receives most of the expelled matter of the disrupted NS. 
Essentially pure BH-torus ejecta of some of our models can reasonably well match the abundance 
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distributions observed in heavy-element deficient metal-poor stars like HD88609 and HD122563 
(right plot in Fig 4). 

Since the major part of the torus ejecta ends up in forming A ^ 130 material, the additional 
mass loss of the merger remnants does not alter event-rate estimates based on comparing yields of 
heavy r-nuclei in the ejecta of the binary-merger phase with the r-process abundances in our Galaxy 
[2, 26]. 
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